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theory of spontaneous emission of a chiral molecule located near a chiral (bi-isotropic) spherical 
particle is developed. It is shown that the structure of photons in the presence of chiral spherical 
particles differs significantly from the structure of TE or TM photons. Exact analytical expres- 
sions for the spontaneous emission radiative decay rate of a chiral molecule with arbitrary electric 
and magnetic dipole moments of transition located near a chiral spherical particle with arbitrary 
parameters are obtained and analyzed in details. Simple asymptotes for the case of a nanoparticle 



^^ \ are obtained. Substantial influence of even small chirality on a dielectric or "left-handed" sphere 

ri" ' is found. It is shown that by using chiral spherical particles it is possible to control effectively the 

cn . radiation of enantiomers of optically active molecules. 
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I. INTRODUCTION 



Chirality is a geometric property of a three-dimensional body not to coincide with its 
reflection in a mirror in any shifts and turns. Such a property, for example, belongs to a 
human hand or a spring. The term "chirality'' was proposed by Lord Kelvin in 1873 to 
explain some special properties of molecules pj. The chirality plays an important role in 
biology and pharmacy, because many complex organic compounds (amino acids, proteins, 
and sugars) have chiral properties. For this reason, a body can react differently with different 
enantiomers of the same substance. For example, the same drug, depending on what type of 
molecules it contains, may have a different taste and smell, or effect differently. In physics, 
there is a considerable interest to chiral media |2|-|5|, in which there is a difference between 
the distributions of left- and right-handed circularly polarized electromagnetic waves. The 
simplest example of such an optically active medium is an aqueous solution of sugar. At the 
present time, the study of the chiral properties is again under special attention because of 
the possibility to create metamaterials based on chiral objects [6|4l2|. 

The spontaneous emission of atoms and molecules in the vicinity of material bodies 

is studied in a large number of works. The influence of dielectric microspheres on the 

spontaneous emission of atoms is considered in 13l-ll6|. spontaneous emission of an atom 

ocated near a microsphere made of material with negative refraction was considered in 

171 ] for first time. Optical properties of chiral spherical particles have also been studied 

quite extensively. At the present time, the scattering of plane electromagnetic waves on 



homogeneous 
microspheres pJ 



and inhomogeneous 



chiral microspheres, including multilayered chiral 



20| . is already considered. The scattering of a Hermitian laser beam on 



a chiral microsphere is studied in [21], the scattering of a plane electromagnetic wave by a 
chiral sphere located in a chiral medium is considered in 22|. The radiation pressure force 
acting on a chiral spherical particle placed in circularly polarized wave is examined in 23 1. 
Finally, the radiation of a point dipole source located inside a chiral spherical particle was 
regarded in [24]. At the same time, as far as we know, there is no investigation devoted to the 
influence of chiral spherical particles on the radiation of chiral (optically active) molecules. 
However, such processes become more and more important because they always appear when 
different biomolecules and drug molecules are investigated and modified by optical methods. 
So, in this paper we will present the detailed investigation of this problem 



The goal of the present work is to develop an approach that allows describing quantita- 
tively the influence of a chiral spherical particle made of an arbitrary metamaterial on the 
spontaneous emission of a chiral molecule. All results will be obtained in the framework of 
nonrelativistic quantum electrodynamics, in the assumption that the radiative linewidth of 
the molecule is much smaller than the frequency of radiation. The last assumption allows 
us to use the perturbation theory. The rest part of the paper is organized as follows. In 
the Section [TTl the canonical quantization of the electromagnetic field in the presence of a 
chiral spherical particle with arbitrary permittivity and permeability is performed. In the 
Section IIIIl we obtain general expressions for the spontaneous emission radiative decay rate 
of a chiral molecule located near a chiral sphere. In the Section [TV] the graphical illustration 
of the results obtained and their discussion are presented. The geometry of the problem is 
shown in Fig. [TJ 

II. QUANTIZATION OF THE ELECTROMAGNETIC FIELD IN THE PRES- 
ENCE OF A CHIRAL SPHERICAL PARTICLE 

To solve the problem of quantization of the electromagnetic field in the presence of a 
chiral spherical particle, let us consider a spherical cavity of an infinite radius A — ;■ oo with 
a perfectly conducting wall and with the a chiral particle located in its center. For simplicity. 



we assume that the nanoparticle is surrounded by vacuum (see Fig.^. T 



le usual procedure 



of quantization of the electromagnetic field in spherical geometry 25l-l27l| is inapplicable in 
the case of chiral particles, and below we will develop a new approach that is vahd in this 
case. 



In the description of classical 
use the method proposed in 



fields in the presence of chiral spherical particles, we will 



28| . At the same tirne, to describe the chiral medium, we 



use the constitutive equations in the Fedorov's form 29| : 



D = £ (E + r/rotE) , 

B = /i (H + r/rotH) , (1) 

where D, E and B, H are the inductions and the strengths of electric and magnetic fields, 
correspondingly, e, n are the dielectric permittivity and the magnetic permeability of the 



chiral medium, and t] is the dimensional parameter of chirality. Substituting equations ([T]) 
into Maxwell's equations, one can obtain 



rotE — iyUxrotH = ikofiH, 

rotH + iexrotE = -ikosE, (2) 

where fco = uj/c is the wave number in vacuum, and x = kot] is the dimensionless parameter 
of chirality. In Eq. ([2]) and everywhere further, the time dependence exp {—iut) is assumed. 
If we rewrite the system ([2]) in the matrix form: 



K = ^-^^ I ^■^'' '" I . (3) 





1 - x^£/U 



-te X£:/i 



then, with the help of the Bohren's transformation for the electric and magnetic fields 



M=J«M. aJ' -"/v^I, (4) 

H y \qr) \ -1,/^i^i 1 

the matrix K can be diagonalized: 

A-'KA = r:" \ I , (5) 

where 




i^L = ^, Kr = — ^, (6) 



are the wave numbers of the left- (L) and right-polarized (i?) waves. The components of the 



transformed field (jl]) satisfy the equations 



28|: 



rotQi = +kLClL, divQi = 0, 

TotQu = -kjiQji, divQ/j = 0. (7) 

Thus, the fundamental solutions of the Maxwell's equations in a chiral medium are the 
left-polarized and right-polarized waves with different wave numbers. 



In the spherical geometry under consideration, in the case of an arbitrary electromagnetic 
field, the expressions for Q^ and Qr can be written down as follows (we will omit a multiplier 
exp {—iut) here and below): 



oo n 



n=l m=—n 
oo n 

Q^ = E E ^^^^ i^^^^^ - ^V'LII) , (8) 



n=l m=—n 



where Amn and Amn are the coefficients of expansions. These coefficients can be found with 
the help of additional conditions. Explicit expressions for spherical vector harmonics N^l>rnn 
and Mil)mn (where J = R, L) are given in the Appendix. 

From Eq. (|H]), it follows that Q^ and Q/j are expressed in terms of fixed combinations 
of vector spherical harmonics Ntpmn + Mi/^mn and Nxpmn — Mif^mn , which always have a 
nonzero component along the radius (see the Appendix) and therefore never can be reduced 
to the usual TE and TM fields. 

In accordance with the general rules of quantization of the electromagnetic field, its ex- 
pansion in the complete system of eigenfunctions of the classical problem (standing spherical 
waves) can be represented as follows (the indices are omitted): 



-^ a^e (s 


,r)-ate^ 
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>r) 
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r) + ath* 


is, 


r) 



E(r) = z^ 

H(r) = ^ "-"^^'^^^"-" ''''\ (9) 

s * 

where a^ and aj are the operators of annihilation and creation of a photon in s-th mode 
respectively, with usual rules of the commutation, e {s, r) and h (s, r) are the modes of 
the electric and the magnetic field, the vector index s = (n, m, v) denotes a chosen set 
of quantum numbers: the orbital (n), the azimuthal (m), and the radial (z/), the asterisk 
denotes the operation of the complex conjugation. In accordance with Eqs. dl]) and ([8]), the 
fields inside the chiral sphere can be written as 



e^" is, r) = Al^l [Nil^i^l + Mt/.^) _ z-^A^^j (iVVL^j - Mi^^^il) 



h'" {s, r) 



/i 



AL^) (AT^^W + M^l^^) - zAl^Ji {N,l^i:>, - M^i^:i) 



(10) 



Outside of the chiral particle, all waves have the same propagation velocity and can be 
represented as 



e°"* (s, r) 



C^NC^^^ +I)(i)MC(') +C(2)Arc(2) +D(2)^^ 

mn >>mn ' van >>mn ' mn ^mn ' mn ^\ 



(2) 



h°"* (., r) = - (C« MC« + D^NCl^ + Cl^lM^l + D^±NC^±) 



:il) 



where Cmn-, Cmn and Dmn, Drrin are some coefficients, and expressions for the spherical 
vector harmonics NQin and MQnn (where j = 1, 2) are given in the Appendix. 

To find the coefficients in Eq. flTTj) . one should use the condition of continuity of tangential 
components of the electric and magnetic fields on the surface of a chiral spherical particle 
having the radius a. Since the molecule is situated outside the particle, we do not need 
to know explicit values of the coefficients Amn and Amn (see Eq. (|H])). For the coefficients 
that describe fields outside of the sphere, one can obtain the following relations from the 
boundary conditions: 






(12) 



in which 



a. 



1-2 



l^n 



An{L)Vn{R)+An{R)Vn{L) 

W„ (L) K [R) + Wn [R] K [L] 

Wn{L) An{R) -Wn{R) An{L) 

Wn (L) v; {R) + iy„ {R) K (L) 

Py„ (L) g„ (i?) + Wn [R] Bn jL) 

W^ (L) Vn (R) + W„ (R) K (L) 
In Eq. ( 1T3|) . we have used the abbreviations (J = i^, -R): 



In 



1-2 



(13) 



iy„ (J) = P^„ (kja) Cl^)' (koa) - ^; (fcja) C^') (koa) , 
K (J) = V'n (A:ja) C^')' (fcoa) - Pi^'n {kjo) (i'^ {k^a) , 
An (J) = Ptpn {kjo) tp'n (koa) - ip'^ {kjo) tpn (koa) , 
Bn (J) = V'n (kja) tp'^ (koa) - Ptp'^ {kjo) ipn (koa) , (14) 

where P = ^/Fjl/fi, ipn (x) and Q (x) are the Riccati-Bessel functions (see the Appendix), 
and the prime near the function means the derivative with respect to its argument. 

To study the interaction between a two-level molecule and the continuum of electromag- 
netic modes modified by the presence of a chiral spherical particle, it is necessary to know 
the density of final states. To find it, one should use the condition of disappearance of 
tangential components of electric field on the inner surface of the cavity. As a result we have 

C« Ci^)' (koA) + Ci^) Cf ' (A:oA) = 0, 
D^M'^ (koA) + D^^Ci'^ (koA) = 0, (15) 

where (n (x) is the Riccati-Bessel function (see the Appendix), and the prime near the 
function means its derivative. 

Substituting Eq. flT2|) into Eq. flTSl) and using asymptotic expressions for the Hankel 
function at A — )■ oo [30], we obtain 

(a„ - (-1)"+^ exp {-2zkoA)) C^± + fi^B^^^ = 0, 

^nCi^l + hn + (-1)"+' exp {-2zkoA)) B^± = 0. (16) 

The system ( IT6l) has nontrivial solutions only if its determinant equals to zero. As a 
result, we have a quadratic equation, which determines two different relationships between 
the coefficients Cmn and Bmn 

^ = -^ i^n + lnT \/{an + In)' ' ^Pl) • (17) 

The relations (TT7|) define two types of electromagnetic modes (photons) occurring in the 
presence of a chiral spherical particle, which are analogous to TM and TE modes in the 



case of a nonchiral spherical particle. The upper sign (— ) in Eq. (fT7|) corresponds to modes 
which we will call A-tjpe photons, while another sign (+) gives us B-tjpe photons. If the 
chirality parameter tends to zero (x — > 0), then, as it follows from Eqs. f ll3p and (I14p . the 
right side of Eq. flT7|) tends to zero or to infinity, which leads to a decoupling between TE 
and TM modes (vector spherical harmonics NCmL and M(^rnn in Eq. (TTT]) ). and as a result 
leads to the appearance of usual TM or TE polarized photons. 

Thus, in the case of a chiral spherical particle we have two types of polarization [A and 
B), which are different from the TM and TE polarizations. 

Substituting Eq. (fT7|l into Eq. ( TT6|l . we find the next asymptotic expressions (A — > oo): 



V c 






where Ug is the mode eigenfrequency, and c is the speed of light in vacuum, and 



Or 



-^(«n + 7n-v/K+^0^^^)- (19) 

From Eq. f llSp . it follows immediately that the density of final states for both A-type and 
-B-type photons has the same value: 

, ^ du A , . 

a [nojs) TTiic 
and does not dependent on the presence of a spherical particle of a final volume in the cavity, 
which is consistent with the Courant's theorem 31|. 

For the normalization of modes, we will assume that there is only one photon in the 
quantization volume. As a result, we obtain the condition: 

^ [dV{d {s, r) e* {s', r) + b (s, r) h* {s\ r)) = 5,,,^u;„ (21) 

Sit J 

V 

where 6ss' is the Kronecker's delta, and the inductions d and b of photonic modes can be 
expressed via the operators of field strengths e and h with the help of the constitutive 
equations ([1]). To calculate the integral in Eq. ( 12T|) . we note that the main contribution to it 
comes from the region r ~ A outside the spherical particle, where there is no chirality and 

8 



where it is possible to use asymptotic expressions for the Hankel functions. As a result, we 
obtain [n = 1, 2, 3, . . .; m = 0, ±1, ±2, . . ., ±n): 

\(r(^)\ \^ = \(r,{2)\ |2_ f^^s (2n + l)(n-m)! 

II mn)A\ II mn)B\ 2Ac^l + \Onf) n (u + 1) (u + my/ ^ ^ 



Thus, expression ( jTTl) . ( 1121) . and ( 12211 . together with the expression for the density of the 
final states (l20l) . fully describe the quantized electromagnetic field in the presence of a chiral 
spherical particle, and allow to describe the interaction of this field with arbitrary atoms 
and molecules. In the next section, we will apply these expressions to find the radiative 
decay rate of the spontaneous emission of a chiral molecule. 

III. RADIATIVE DECAY RATE OF SPONTANEOUS EMISSION OF A CHIRAL 
MOLECULE LOCATED NEAR A CHIRAL SPHERICAL PARTICLE 

The radiative decay rate of the spontaneous emission (the radiative linewidth 7) of a 
chiral molecule located near a chiral spherical particle can be found by making use of the 



Fermi's "golden rule" 



32|: 



7 = — ^ \{imtial\Hint\ final) f p (u) , (23) 

final 

where summation is performed over all possible final states, which include both types of 
photons {A and B, see Eq. (ITTj) ). p (w) is the density of the final states (l20l) . and Hint is the 
interaction Hamiltonian of a chiral molecules with the electromagnetic field. In this work, 
we assume that the Hamiltonian is Hermitian, i.e. there are no losses in the material of a 
particle. If there are losses, there is a possibility of nonradiative transition from the excited 
state to the ground state. The expression for the nonradiative part of the spontaneous 



decay rate can be obtained within the 
the nonstandard quantization scheme 



ramework of the classical approach 33| or by using 



34 



37|. In further calculations, for brevity we will 



consider only one channel of the decay of the initial state of the molecule (transition e ^ g), 
i.e. the two-level molecule. To take into account the possibility of transition of the molecule 
into several states, it is enough to sum the partial linewidths. 

The Hamiltonian of interaction of a chiral molecule with the electromagnetic field in 
Eq. (|23|) can be written down in the next form ^ 



Hint = - f dE (ro)) - f liiH (fq)) , (24) 



where Fq is the radius vector of the molecule position, d and m are operators of the electric 
and the magnetic dipole moments of the molecule, and field operators E and H are defined 
in Eq. ([9]). We suppose that the initial state corresponds to the vacuum state of the field 
and to the molecule in the excited state, {initial\ = {e\{vac\, and that the final field state 
correspond the single photon (of A-tjpe or B-type) and the molecule in the ground state 
\final) = \ln.m,u)\g)- Besides, for definiteness we will consider a spiral model of a molecule 
in which electrons are constrained to move along a helical path. Substituting Eq. ( 12^ into 
Eq. fl23|) . we obtain the following expression for the total radiative decay rate of spontaneous 
emission of the molecule: 

7eg = 7eg + 7fg, (25) 

where, for example, in the case of photons of A-tjpe 



4 = J-P (^) E I (do< (^' ^' ^' ^o)) - (^o^A * (ri, m, u, Fq)) | ', (26) 



n,m 



where (e 



g) = do and (e|m|(7) = — imo are dipole moments of the considered transition 
of the molecule on the frequency u ^ Us (see Eq. (IT8|) ). The chosen phase difference between 
the electric and magnetic dipoles is due to the fact that the magnetic moment operator is 
purely imaginary. This definition also agrees with the adopted spiral model of a molecule. 
For more complicated chiral molecules, there can be a different phase relation between 
electric and magnetic dipole matrix elements. 

Let us, for clarity, call the molecules with parallel do and mo the "right" molecules, while 
the molecules with antiparallel do and mo will be referred to as the "left" molecules. Explicit 
expressions for the component of Eq. (l25l) for the mode of photons of i?-type one can obtain 
from Eq. f l26|) by replacing the index A — ?■ 5. 

Let us assume for definiteness that the molecule is located on the z-axis of a Cartesian 
coordinate system at the point tq > a. In this case, only components with m = 0, ±1 are 
nonzero, and the explicit expressions for Eq. ( I25l) will have the following form: 



7e1 = it' + <' + 7e1'°, (27) 

10 



where 



<'"' =^E T-^K^o. - ^doy) « (Vo) + T„\(i)' (Vo)) 
-On {doy + ifiox) (^n {koro) + L;^C^^^ (A^o^o)) 

- (moj, + irriQ^) (V^„ (A;oro) + T^C^ {hr^))\\ (28) 



oo 



- {do, + idoy) « (A;oro) + T^C^' {hro)) 

+ {rriQy - irriQ,) (?/'„ {k^ro) + T^Ci^^ (A^o^o)) 

-0„ (mo. + ^mo,) (^; (A;oro) + L^C^^^' {k^ro)) \\ (29) 



and 



7, 



2 °° 



Ao _^^ (2n + l)n(n+l) ^^^^ ^^^ ^^^^^^ ^ ^^^^(,) ^^^^^^^ 



eg 

ii.ri.i\i /-, I I If I 

n 



+0„mo. {^n (koro) + L^Ci'^ {koro)) f . (30) 

In the case of photons mode of 5-type, one can obtain 



75 = 7er'+7f/+7f/, (31) 



where 



B,-l _ '^1 



^"^ 2hr, 



2^^ ^^^r^l<^" (^ox - idoy) « (Vo) + 7;f d^^' (Vo)) 

0„=l ^+ \On\ 



- {doy + ido,) {ipn (koro) + L^Cn^ (/co^o)) 

+ (mox - imoy) « (koVo) + L^Cn^' (koro)) 

-On {moy + tmo,) {^n (koro) + T^C^ (koro)) \\ (32) 



11 



B 1 ko ^ 2n + 1 



+0n {moy - imoxO (V^n (A;o^o) + T^Cn^ (A^o^o)) 

- (mo. + zmoy) (^; (fcoro) + ^f Ci'^' (/^oro)) f, (33) 



and 



7, 



B,0 



2 ^ (2n + l)n(n + l) 



E ^ ,1^2 > A i^n (fcoro) + Ti^e^ (fcoro)) 



+mo. (V^n (koro) + L^C^ (A;oro)) f . (34) 

In Eqs. (128])-([30D, ([32])- ([MD, we have used the following notations: 



T^ = 4 ( «n - 7n + V («" + 7n) ' ^Pl " 2 1 , 
Li = \ [in - «n - V("" + 7n)' " 4/3^ - 2^ , 



T^ = \\an-ln- 7(a„ + 7n)'-4/32 - 2 

^^f = ^ (7™ - ttn + 7(a„ + 7n)'-4/32 - 2") . (35) 



In the case when the chirality parameter of a spherical particle is equal to zero, Eqs. fl28 
(130|) . (!32|) - (!M|) are simplified and have the following form: 



00 



it" =^E (2^ + 1) K^o^ - ^^0^) K (^0^0) + ^'^C"'^' (^0^0)) 



2^ 1 



- (moy + ^moo;) (^n (^0^0) + T^Cn^ (koro)) \ , (36) 



kn 



<' =mhll ^2r^ + 1) \{do. + ^doy) (^; (A:oro) + T^^C^' (^oro)) 

- (moy - imox) {'4'n (A^o^o) + T^Cn^ (koro)) \ , (37) 

12 



|2 oo 

^er = TiAE (2^ + 1) n (n + 1) |V^„ (fcoro) + T^^^ {k,r,)\\ (38) 



Afi ^ ^Mozl 



and 



, oo 

^ n=l 

- (mox - imoj;) (V^^ (fcoro) + L^Cn^' (A^o^o)) | , (39) 



, oo 

^ n=l 

- (mox + imoy) « (fcoro) + Lfxll^' (Vo)) f , (40) 



91 P °° 
lef = 4?AE (2^ + l)n{n + l) |^„ (Vo) + I^^d') ( Vo) r, (41) 



'eg 



^^oro" tt 



where 



rpA 

n 



B 



Pipn [koay/ejl) Cn^^' (koa) - ip'^ [koa^/ejl) Cn^^ (koa) J ' 
tpn {kpay/eji) j)'^ {kpa) - Pip'^ {kpa^/ejl) ipn jkoa) \ 
^tpn (koay/ejl) Cn {kpa) - P'ip'n [kpa^jejl) Qt {kpo) ) 



(42) 



Eqs. (!36|) -( l4T]) correspond to components of the radiative decay rate of spontaneous emis- 
sion of a chiral molecule located near the dielectric sphere. If we put mo = or dp = in 
these expressions, it will be possible to find known relations for the radiative decay rate of 
spontaneous emission of an atom located near a dielectric spherical particle, and which has 
nonzero electric or magnetic dipole moment, respectively 



Note, that despite of the fact, that this derivation was performed in assumption of absence 
of losses in the material of the particle, the expressions obtained will also exactly describe 
radiation losses for nanoparticles made of real materials, i.e. for complex e and yU. It follows 
from the fact that classical and quantum calculations give the same analytical expressions 
for the radiative decay rate of spontaneous emission through the Green's function of a 

13 



classical problem J40|, |4l|. Since the classical approach is valid in the case of both absorbing 
and nonabsorbing bodies, this means that for a generalization of the expressions obtained 
for the case of a real (absorbing) material of a spherical particle, it is only necessary to 
substitute appropriate complex values of dielectric permittivity and magnetic permeability 
into obtained expressions. 

To calculate the relative radiative decay rate of spontaneous emission of chiral molecules 
located near a chiral spherical particle, one should normalize the expression (1251) to the 
value of the radiative decay rate of spontaneous emission of the molecule in the absence 
of the particle. For this purpose, we use Eqs. (15^ -( HT]) where we put T!^ = L^ = 0. 
By substituting the expressions obtained in such a way into Eqs. (127|) . (13T|) and (125|) . and 
performing summation, it is possible to find the radiative decay rate of spontaneous emission 
of a chiral molecules in free space: 



70,. ^ f IdolV f Kll (43) 



This expression, of course, coincides with the sum of the known expressions for the linewidths 
of electric and magnetic dipoles 



39[ |. because in the absence of chirality the interference of 



radiation does not occur. 



IV. ANALYSIS OF THE RESULTS OBTAINED AND ILLUSTRATIONS 

Cumbersome formulas obtained in the previous section cannot fully reveal the features 
of the spontaneous decay of a chiral molecule in the system under consideration. In this 
section, we will explore the expressions found in the typical special cases. 

In a very important case of a chiral nanoparticle, i.e. a particle having a size much smaller 
than the wavelength of the radiation, the main contribution to Eqs. (I28l)-(l30l). (!32l) -(1Ml) is 
from term with n = 1. Expanding the coefficients of Eq. ( IT3|) in a series of koa — )■ and 
keeping only the principal terms, we obtain 

14 






3 ^ " ^ V(^ + 2)(/x + 2)-4x%/x 
Pi ~o^ (^oo) 



The condition that the denominator of Eq. (1441) will be a zero 



{e + 2) (/i + 2) - Ax^e^ = (45) 

corresponds to the condition of a chiral-plasmon resonance occurring in a spherical nanopar- 
ticle. The roots of Eq. (H5|) determine values of the constants e and /x, that correspond to a 
chiral plasmon. It should be noted, that the resonant values of e and /i are not independent, 
and their values differ from the corresponding values in the case of plasmon resonance in a 
nonchiral particle, where ( H5|) is split into two independent equations, e + 2 = and /i + 2 = 



42| 



To find an expression for the spontaneous emission radiative decay rate in the case of a 
nanoparticle, one should expand n = 1 terms in Eqs. (l28l) -( !30l) . ( !32|) -(!34l) into a series of 
k^a — )• and keep only the principal terms. However, it is a too long procedure, and here 



we will apply a different approach used in 43|. If the nanoparticle and the molecule are 



placed close enough, then their total radiation will be determined by their total electric and 
magnetic dipole moments: 

7e. = ^|do + 5df + ^l-^mo + 5m\\ (46) 

where Sd and (5m are the electric and the magnetic dipole moments induced in the nanopar- 
ticle while do and — imo are the oscillating dipole moments of the molecule. Explicitly, these 



values can be written as follows |43J : 



5d = ttEE^o (ro) + "shHo (ro) , 
(5m = ohe'Eo (ro) + a^nHo (ro) , (47) 



where the polarizability of a chiral spherical nanoparticle in a homogeneous field is 
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3(£-l)(/i + 2) + 2xV 



OlHH 



3(^- 


2)(/x + 2)-4x%/i 
- 1) (e + 2) + 2x V 


3 


2)(/x + 2)- 


- 4x^£:/i 



a^ff = -aHS = a - — — -- — — -—, — . (48) 

[e + 2) (/i + 2) - 4x%/i 

The fields in Eq. ( I47l) have the form: 

3n(ndo)-do .3n (nmo) - mo 

Eo(r) = ^^ , Ho(r) = -z -^ , (49) 

where n is the unit vector in the direction from the center of the particle to the observation 
point. 

In practice, the orientation of the molecule can be arbitrary with respect to the nanopar- 
ticle surface, therefore to get an effective radiative decay rate of spontaneous emission one 
should average Eq. (H6|) on the orientations of molecules, or, equivalently, on the unit vector 
n. As a result we obtain 

Ak^cf (2 2 1 

611 I rQ Tq J 

where ^ is defined by the relation mo = .^do. 

As a rule, the magnetic dipole moment of the molecule is much smaller than the electric 

dipole moment |mo| ^ |do|- The chirality parameter, even in hypothetical metamaterials, 

is also small {x ^ 1)- This fact determines that the second term in Eq. f l50|) . corresponding 

to the induced electric dipole moment, is usually greater than the term corresponding to the 

induced magnetic dipole moment. Thus, the effective interference between the electric and 

magnetic fields is possible only if the following two conditions take place: 

1. In the system under consideration, a chiral-plasmon resonance must be present, i.e. 
the condition (^^ must be satisfied. Under this condition, the contribution of a magnetic 
radiation increases. 

2. The electric dipole moment induced in the nanoparticle should be zero, i.e. the next 
condition must be satisfied (see Eq. ( fJOjj j 

OiEE - iioiEH = 0. (51) 
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The solution of the system of Eqs. (H5|) and (15T|) determines the values of dielectric 
permittivity and magnetic permeability of the nanoparticle, which correspond to minimal 
values of the radiative decay rate of a chiral molecule. It means that the interference between 
the electric dipole and magnetic dipole radiation becomes maximal and destructive if 



^l^ 



2rfn 



do + 2xmo ' 



e -)■ -. (52) 

niQ + 2xdo 

From Eq. ( 152|) . it follows that by changing the sign in ttiq, i.e. when the chirality of the 
molecule is changed, the "resonant" magnetic permeability jl varies only slightly (x ^ 1, 
mo ^ do) and approximately equals to —2. On the other hand, the "resonant" dielec- 
tric permittivity e may have different signs for molecules with different chirality. This 
means that both nanoparticles with simultaneously negative e and yU (double negative - 
DNG-metamaterials), and nanoparticles with negative /i and positive e {fi negative - MNG- 
metamaterials or magnetic plasma) are suitable for the effective control of radiation of chiral 
molecules. From the viewpoint of practical implementation of such nanoparticles, the most 
suitable are nanoparticles with a negative /x and positive e, because they can be implemented 



by using a well developed technology of Split Rings Resonators [44 |45|. On the other, the 



chiral DNG-metamaterials can be also synthesized |12 |. 

The above expressions give a good description of properties of arbitrary chiral molecules 
radiation of near nanoparticles, i.e., when the retardation can be ignored. In the case when 
retardation effects are significant, it is necessary to use the full expression for the decay rates 
given in the previous section. 

In any case, the process of spontaneous decay of a chiral molecule located near a chiral 
spherical particle is very complex and below for clarity we will give graphical illustrations 
of some possible regimes of interaction of a chiral molecule with nanoparticles of different 
composition but neglecting losses in them. 

The Figure [2] shows the distribution of absolute values of the coefficient T/ (see Eq. ( I35l) ) 
depending on the dielectric permittivity and the magnetic permeability of the material 
of non chiral (Fig. |2]^a)) and chiral (Fig. Mjo)) spherical particles for different values of koa. 
Large values of this coefficient correspond to a resonance. This figure clearly shows that there 
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is a significant difference between these cases, as refiected in the fact that nonzero chirahty 
of the particle leads both to an increase in the number of observed modes (whispering gallery 
(WG) modes in dielectric "right-handed" (RH) material and in "left-handed" (LH) material, 
LH surfaces modes [171]), and to the change of their structure as compared with a nonchiral 
particle. This fact is due to the simultaneous presence of both left- and right-polarized 
waves in a chiral particle, which in turn leads to excitation of resonances analogous to TM 
and TE resonances of nonchiral sphere, simultaneously. The Figure [2] is very important 
because it allows one to qualitative understand the structure of the radiative decay rates of 
spontaneous emission of a chiral molecule placed near the different nanoparticles. 

It is very interesting that even for a dielectric sphere with a small admixture of chirality, 
strong and unexpected effects occur. In Fig. [3l the radiative decay rate of spontaneous 
emission of a chiral molecule placed near a chiral spherical particle with positive values of 
permittivity and permeability is shown as a function of the size parameter. As it is seen in 
this figure, the presence of chirality in a spherical particle leads to a shift of the maxima of 
the spontaneous emission radiative decay rate in comparison with the case of a nonchiral 
particle. Indeed, as it is known, a condition that allows to estimate the position of the 
resonance is given by koay/ejl ^ n + 1/2, where n is the orbital quantum number. In the 
case of chiral spherical particles, the left- and the right-polarized waves exist simultaneously, 
and kia > kQa^/eJl if x > 0. Hence the above condition should be changed to k^a ~ n + 1/2. 
As a result, this leads to smaller resonant values of k^a than for a nonchiral particle (see 
Fig. [2]). Even more interesting feature of ordinary dielectric particles with a small admixture 
of chirality is a substantial increase in the quality factor of whispering gallery modes. This 
figure clearly shows that the corresponding linewidth can be increased by the factor 9 or 
even more. 

Note also that in Fig. [3], only a slight difference between the radiative decay rates of 
spontaneous emission corresponding to the different orientations of the transition magnetic 
dipole moment of a chiral molecule at a fixed electric dipole moment is observed. However, 
one should expect a growth of this difference in the case of a gradual increasing of absolute 
values of the magnetic dipole momentum. 

Figure H] shows the radiative decay rate of spontaneous emission of a chiral molecule 
located near chiral spherical nanoparticles [k^a = 1) as a function of the chirality x = ^ol- 
As it is clearly seen in this figure, changing the chirality of the sphere has the greatest impact 
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on the rate of spontaneous decay of the molecules near dielectric and left-handed spherical 
particles, in which high-Q modes can be excited. When the chirality parameter approaches 
to the critical value Xcru = '^1 \/^i then the number of oscillations in the radiative decay 
rate increases rapidly, and fc^ — )■ oo (see Eq. ([6])). In the case of metallic particles, the 
dependence of the radiative decay rate of spontaneous emission on the parameter x, which 
is shown in Fig. HJ is weak, due to imaginary values of ^/eJl, and due to the absence of 
propagating waves. 

Chiral molecules play an especially important role in biology and pharmacy. Therefore, 
it is extremely important to arrange the separation of "right" and "left" enantiomers of 
molecules in the racemic mixtures. This can be done in various ways, for example by 
using radiation pressure forces in the electromagnetic field of left- or right-handed circularly 



polarized electromagnetic waves 23|], or by using spiral optical beams 46|. In addition, such 
selection can be efficiently performed with the help of chiral nanoparticles 43 1 . 

Figure [5] shows the radiative decay rate of spontaneous emission of a chiral molecule 
located near the LH chiral spherical nanoparticles as a function of the dielectric constant 
for a fixed value of magnetic permeability. As it well seen in the figure, for molecules that 
differ only in the orientation of the magnetic dipole moment of the transition (the "left" 
molecule with mo = — O.ldo, and the "right" one with hiq = +0.1do), the radiative decay 
rate takes significantly different values, and thus, under such conditions it is possible to 
control effectively the radiation properties of enantiomers. It is important to note that both 
the results of calculations obtained in the framework of the full QED theory (present paper), 
and in the framework of a simplified quasistatic calculations presented at the beginning of 



this section (see also 42]) yield almost identical results, confirming the correctness of both 
approaches. 

The Figure [6] shows the ratio of the effective radiative decay rate of spontaneous emission 
of the "left" molecule and the radiative decay rate of the "right" molecule and vice versa. 
From this figure, it follows that if the condition f l52|) is satisfied, then the decay rates of "left" 
and "right" molecules differ by the factor 15 or 60 or even more, depending on the chirality of 
the molecule considered as a reference one. In other words, nanoparticles with the parameter 
given by Eq. ( I52l) will enhance the radiation of the "right" molecules and slow down the 
radiation of the "left" molecules and vice versa. Let us stress that "left-handedness" of a 
chiral sphere or its negative /x are of crucial importance for such discrimination. Possible 
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applications of the effect of discrimination of the radiation one can find in 43 1. 



V. CONCLUSION 

Thus, in the present work the analytical expressions for the radiative decay rate of spon- 
taneous emission of a chiral molecule located near an arbitrary spherical particle with chiral 
properties were obtained and investigated within full QED theory. Using this approach, the 
spontaneous emission decay rates of chiral molecule placed near a spherical particle made of 
different materials (dielectrics, metals, "left-handed" metamaterials etc.) were investigated 
in details. It is shown that nonzero chirality of a spherical particle leads to an increase in the 
number of excited modes in comparison with nonchiral one due to coupling between left- and 
right-polarized waves. Quasistatic expressions for the radiative decay rate of spontaneous 
emission of chiral molecule located near a chiral nanosphere were obtained and investigated. 
These expressions are in good agreement with the exact QED results and allow one to es- 
timate the parameters of nanoparticles when the radiation of "right" or "left" molecules is 
suppressed. It was found that for the suppression of the radiation of the "left" molecules 
the sphere should be made of MNG metamaterials, while for the suppression of radiation 
the "right" molecules the DNG materials should be used. 

The results obtained can be used to calculate the radiative decay rate of spontaneous 
emission of chiral molecules in the vicinity of the chiral spherical particles, for interpretation 
of experimental data on interaction of chiral molecules and particles, and for detection and 
selection of chiral molecules with the help of chiral particles. 

Appendix: Spherical vector harmonics 

Spherical vector harmonics that describe the electromagnetic field inside a chiral spherical 
particle have the form: 



^. ijl,{kjr)mP;^ {cos 6) ^,^^^ 
kjr sm & 
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(A.2) 



where O<r<oo,O<0<7r and < < 27r are spherical coordinates, e^, eg, e^ are 
unit vectors of the spherical coordinate system, ipn {kjr) = {TTkjr/2) ' Jn+1/2 (kjr) is the 
Riccati-Bessel function |301], where Jn+1/2 (kjr) is the Bessel function [30|, the prime near 
the function means the derivative with respect to its argument, P™ (cos 6*) is the associated 
Legendre function 30|, and the index J = L, R. 

Spherical vector harmonics that describe the electromagnetic field outside a chiral spher- 
ical particle have the following form: 
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(A.3) 



'(i) 



^(i) 



k^r sva.0 k^r oO 



where (n (kor) = {Trk^r /2) ' H^li/2 (^o^) is the Riccati-Bessel function 30|, where -f^„+i/2 (^o^) 



(A.4) 



is the Hankel function of the first (j = 1) and second (j = 2) kind, respectively, and the 
prime near the function means the derivative with respect to its argument. 

More information about the properties of spherical vector harmonics can be found, for 
example, in J47|. 
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FIG. 1. The geometry of the problem of spontaneous emission of chiral molecules located near a 
chiral spherical particle. 
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FIG. 2. (Color online) The dependence of the absolute values of the coefficient T4 on the permit- 
tivity (e) and the permeability (fi) of chiral spherical particle for different values of koa. (a) The 
chirality parameter x = 0- (b) The chirality parameter x = O-l- The particle is placed in vacuum. 
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FIG. 3. (Color online) Radiative decay rate of the spontaneous emission of a chiral molecule located 
in close vicinity to the surface of a chiral dielectric spherical particle (rg — )• a) with e = 6 and fi = 1 
as a function of its size, koa. (a) The transition electric dipole moment of the molecule is oriented 
along the x-axis (dg = dQBx), i.e. it is tangential to the surface of the particle, (b) The transition 
electric dipole moment of the molecule is oriented along the z-axis (dg = doez), i.e. it is normal to 
the surface of the particle. The solid line corresponds to the transition magnetic dipole moment of 
a molecule oriented along the x-axis, and the dashed line corresponds to the transition magnetic 
dipole moment of a molecule oriented along the z-axis. The value of the transition magnetic dipole 
moment {mox,iT^Qz} = O.ldo- The particle is placed in vacuum. 
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FIG. 4. (Color online) The radiative decay rate of the spontaneous emission of a chiral molecule 
located in close vicinity to the surface of a chiral spherical particle (rg — )• a) with the size koa = 1, 
as a function of the chirality parameter x- A dielectric sphere has e = 4 and /x = 1, a metal sphere 
has e = —4 and ;U = 1, an LH sphere has £ = —4 and ^ = —1.11. (a) The transition electric dipole 
moment of the molecule is oriented along the x-axis (dg = doex), i.e. it is tangential to the surface 
of the particle, (b) The transition electric dipole moment of the molecule is oriented along the 
z-axis (do = doBz), i.e. it is normal to the surface of the particle. The solid line corresponds to the 
transition magnetic dipole moment of the molecule oriented along the x-axis, and the dashed line 
corresponds to the transition magnetic dipole moment of the molecule oriented along the z-axis. 
The value of the transition magnetic dipole moment {m-oaiW-Oz} = 0.1 do- The particle is placed in 
vacuum. 
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FIG. 5. (Color online) The radiative decay rate of the spontaneous emission of a chiral molecule 
located in close vicinity to the surface of a chiral LH spherical nanoparticle (ro — )• a) with e = 
e' + iO.l, II = —1.6, X = 0.2, and kQa = 0.1. (a) The electric and magnetic dipole moments of the 
transition of the molecule are oriented along the x-axis (tangentially to the particle surface), (b) 
The electric and magnetic dipole moments of the transition of the molecule are oriented along the 
z-axis (normally to the particle surface). The dashed line shows the asymptotic expression (146 p . 
The nanoparticle is placed in vacuum. 



28 



fJ^(Mt)/jf^ihght) 




yf^ihghtyfJ^iMt) 




FIG. 6. (Color online) (a) The ratio of the effective radiative decay rate of spontaneous emission of 
"left" molecules with niQ = — O.ltio to the effective radiative decay rate of spontaneous emission of 
"right" molecules with niQ = +0.1do and (b) vice versa as a function of the dielectric permittivity 
(e) and the magnetic permeability (/i) of the material from which chiral spherical nanoparticle is 
made. The chirality parameter x = 0.2. The molecule is placed in close vicinity to the surface of 
the spherical nanoparticle (ro — )• a). The white line shows the position of chiral-plasmon resonance 
in the nanoparticle (see Eq. (|45p ). The nanoparticle is placed in vacuum. 
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